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Abstract. We examine the time dependent amplitude (j)j (t) at each vertex j of a 
continuous-time quantum walk on the cycle Cn- In many cases the Lissajous curve of 
the real vs. imaginary parts of each tpj {t) reveals interesting shapes of the space of 
time-accessible amplitudes. We find two invariants of continuous-time quantum walks. 
First, considering the rate at which each amplitude evolves in time we find the quan- 
d(t)j (t) 
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is time invariant. The value of T for any initial state can be 
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minimized with respect to a global phase factor to some value An operator 

for Tmin is defined. For any simply connected graph g the highest possible value of 
Tmin with respect to the initial state is found to be T™“f = () where Xmax is the 
maximum eigenvalue in the Laplace spectrum of g. A second invariant is found in the 
time-dependent probability distribution Pj (t) = \4>j (t)| of any initial state satisfying 
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of n vertices. 
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1. Introduction 

Quantum walks on graphs come in two different versions, continuous-time and discrete¬ 
time. For a comprehensive review of quantum walks see Venegas-Andraca [1]. The 
main interests in quantum walks lie with the development of quantum search algorithms 
[21 |3l II] . In this paper, we concentrate on the amplitudes of a continuous-time quantum 
walk (CTQW) on connected graphs. A graph G = (V, E) is defined by the set of vertices 
V which are connected by the set of edges E. The adjacency matrix A and degree matrix 
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D which describe the graph are defined as follows: 

1 if there is an edge connecting vertices j and k 
0 otherwise 

deg {vj) ii i = k 
0 otherwise 

where deg (vj) is the degree of vertex Vj. 

Generally, the nnitary CTQW operator U{t) is defined in terms of the Laplacian 
matrix L of the graph 015] 




(la) 

(lb) 


L = D — A (2a) 

U it) = (2b) 

For k-regnlar graphs, the operator is conventionally defined in terms of the normalized 
adjacency matrix alone [f] [Zl IH| 

U {t) = (3) 

The qnantnm walker is described by a time-dependent qnantnm state |0 (t)). The initial 
state vector |0(O)) is an ordered set of components, each component corresponding to 
the initial amplitude at each vertex of the graph (we will label each vertex starting with 
0 ) 


10 (0)) = [00 (0), 01 (0), 02 (0),..., 0„_1 (0)] (4) 

such that 

n—1 

(<i'(0)l</.(0)>=5^fe(0)|" = l. (5) 

3=0 

The time development of the quantum walk then becomes 

f/(f) 10(0)) = 10(f)), (6) 

and the time-dependent amplitude at vertex j is 

0jW = O10W)- (7) 

2. The time-dependent amplitudes at each vertex of Cn 

A cycle graph Cn is a 2-regular graph with n vertices. The corresponding U (t) will be 
an n X n circulant matrix completely defined by an n component circulant vector a («)0 

U (t) = CIRCn [tto(^) (t)^a2 (0 , ■ ■ ■ .On-l (f)] - (8) 


f This essentially amounts to a change in time-scale and a time-dependent global phase factor. 
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For any cycle Cn the components aj {t) are given by Ben-Avraham, et al. [8] as 
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^ n—1 

(t) = i . (9) 

^ k=0 

This is eqnivalent to the definition in eqnation The amplitndes aj (t) are periodic 
only for n G {1,2,3,4, 6}. This is dne to the fact that all the cos (27r|) terms in 
eqnation are rational only when n G (1, 2, 3,4, 6}. Thns, is the smallest cycle with 
aperiodic amplitndes. When the qnantnm walker is initially localized at vertex j = 0, 
i.e., 00 (0) = 1, the time-evolved amplitnde at each vertex is eqnal to the corresponding 
component of the circnlant vector, 0^ [t) = aj {t). A plot of the real and imaginary parts 
of 00 (t) is given in Fig. There is no obvions phase relationship between the real and 
imaginary parts; however, the conventional tool for investigating phase relations between 
two different waveforms is the Lissajons cnrve or parametric plot. The Lissajons cnrve 

Re[0o(t)] , lm[0o(t)] 



Figure 1. A plot of the time-dependent real and imaginary parts of 0o (0 on a cycle 
Cs. 

of the real vs. imaginary parts of each amplitnde is presented in Fig. Becanse the 
0j (t) are aperiodic, as f —)■ cx), the Lissajons cnrves will become densely filled-in and 
the shape defines a continnous space of accessible amplitndes. A plot of the probability 
at vertex j = 0, Pq (t) = |0o (t) P is given in Fig. The greatest probability in the 
time interval 0 < t < 100 is Pq ~ 0.996259 at t ~ 72.9715. 

When n G (even integers} eqnation yields amplitndes which are pnrely real 
or pnrely imaginary. In these cases, to create Lissajons cnrves with descriptive phase 
relationships we introdnce a time-dependent global phase factor. There is a snbjective 
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n = 5, j = 0 n = 5, j = 1 n = 5, j = 2 n = 5,j = 3 n = 5,j = 4 

Im Im Im Im Im 



Figure 2. The Lissajous curves of the time-dependent amplitude at each vertex 
of Cs when the quantum walker is initially localized at vertex j = 0, plotted from 
0 < t < 200. (The unit circles are included for scale. The green and red dots indicate 
the initial and final values in each Lissajous curve.) 


I0o(t)|" 



Figure 3. The plot of j^o (t) P corresponding to Fig. from 0 < t < 100. 
100 (t) P 0.996259 at t ~ 72.9715. 


component to what the “right” phase factor should be. By trial and error we hnd that 
a global phase factor of a{t) = has an evocative effect. The Lissajous curves for 
Cg with initial state 0o (0) = 1 is presented in Fig. Plots of |0o (t) P and |04 (t) p are 
given in Fig. 

Inspection of Figs. and suggest that the quantum walker will nearly return to 
its initial state at an aperiodic rate. Likewise, inspection of Figs. and suggests the 
quantum walker will essentially oscillate (though aperiodically) between vertices j = 0 
and j = 4. It is an open question if there exists a time t > 0 such that |0o (^) P = 1 
when 100 (0) P = 1 for any Cn when n 0 {1, 2, 3,4, 6}. The problem of recurrence of the 
initial state of the quantum walker in a discrete quantum walk on Cn has been solved 
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n = 8, j = 3 



n = 8,j = 7 



Figure 4. The Lissajous curves of a CTQW on Cg. The quantum walker is initially 
localized at vertex j = 0. Note that the amplitudes at vertices j = 2,6 are periodic. 
(Equationis augmented with a time-dependent global phase factor of 


\<l>o{t)d, I04(t)|^ 



Figure 5. Plots of \(j)o (t) P and \4>a (t) P corresponding to Fig. [^from 100 <t< 200. 
\(j)o (t) P ~ 0.995731 at t ~ 106.722 and \(j )4 (t) p ~ 0.999633 at t ~ 182.185. 
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by Dukes [lO]. 

We have inspected the amplitudes of a CTQW on through n = 201. With 
increasing n the amplitudes become distributed over a greater number of vertices. 
Consequently, the area covered by the Lissajous curves become progressively smaller 
as n increases, also the shapes are often surprising. 


3. A kinematic invariant of a CTQW 


The rate at which the quantum walker state vector |0(p) changes in time is 

i4)> = 

= jU (t) |^( 0 )> 


( 10 ) 


Because an operator commutes with a function of itself and L is real-Hermitian we arrive 
at 

n—1 

(0(t)l0(t)) = 

j=0 

= (0(0) I ^^10(0)) 

= T. 

It is tempting to think of T as the “total kinetic energy” of the CTQW on a graph. 
The value of T is not invariant to a time-dependent global phase factor. Taking 
|0'(^)) = e*®*|0(t)) we readily obtain 


d(t)i it) 
dt 


(4|0'(,)) = (0(o)|(L-0)V(o)) 

= T{e). 


( 12 ) 


The Hermitian matrix L has a full set of orthonormal eigenvectors 
{|uo), Id), •••, |W-i)} with corresponding eigenvalues {Aq, Ai, ..., A„_i}. The initial state 
10(0)) can be expressed as a superposition of the eigenvectors 

10(0)) = ®o|w) + flll^^l) + ■ • ■ + Un-llw-l)- (13) 

In this representation T (6) becomes 

T (9) = (0(o)|0(o))6*^ — 2(0(0) |T|0(o))6' -1- (0(o)|T^|0(o)) 

( n-l \ /n-1 \ /n-1 \ 

j 6^^ - 2 I j ^ + [ ) (14) 

j=0 ) \j=0 ) \j=0 / 

= + B9 + C. 


T ( 6 ) is a convex parabola with a vertical axis of symmetry. Plots of T ( 6 ) for a few 
selected initial states on C 4 are shown in Fig. For any initial state there will be a 
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J{B) 



Figure 6. Plots of T [ 6 ) for 6 different initial states on the cycle graph C 4 . 


minimum value T^in = T {6y) where 6 y is the abscissa value of the coordinates for the 
vertex of the parabola. A fundamental relation in geometry gives the coordinates of the 


vertex in terms of the coefficients of equation 14 


_ ( B 52 - 4AC'\ 

“ V id } 

= ((0(0) I-^10(0)), (0(0) I-^^10(0)) - (0(O)|-^|0(O))^) 


(15) 


where we have A = (0(o)|0(o)) = 1- 

We hnd that T^m can be represented as an operator Tmin acting on the idempotent 
density matrix p representing the initial state 


Tmin p = L^p- LpLp 
= L^pp — LpLp 
= L{Lp- pL) p 
= L [L, p] p, 


(16) 


so that Train = Tr (L [L, p] p). 


4. The extrema of Tmin with respect to the initial state 

Train IS a measure of the total “dynamism” of the amplitude \4>{t)) minimized with respect 
to a global phase factor In other words, the dynamism is due to the time dependent 
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norm of each component of \4>{t))- Thus, the stationary states of the operator U{t) will 
have the lowest possible value = 0. These states correspond to the normalized 
eigenvectors of the Laplacian matrix L. A stationary state, as is indicative of the name, 
will have a constant, time-independent, norm and thus the probability distributions over 
the vertices of the graph will be constant in time. Initial states with values of T^m > 0 
will produce probability distributions which vary in time. Higher values of Tmin are due 
to the range in norm values as well as their rate of change. 

The maximization of T^m, Tmtny does not define a unique initial state. The 
eigenvalues of the Laplacian matrix L can be put in increasing order 


0 — Ao < Ai < A 2 < ■ ■ ■ < An-i. 


(17) 


If the graph is simply connected, Ai > 0. If the greatest eigenvalue A„_i = Xmax is 
degenerate with multiplicity m then 


Xmax 


X 


n—1 


X 


n—2 — 


^n—rri' 


(18) 


The initial states which maximize Tmin are a superposition of the eigenvectors associated 
with the eigenvalue 0 and an arbitrary unit vector in the eigenspace of Xmax of the form 


Then 


10(0)) = {ai\Vn-l) + a2\Vn-2) H - h am\Vn-m)) , 

|q;iP -|- |Q(2p + ■ ■ ■ + |Q(mP = 1- 


rj-^r 


(0(O)|T^|0(O)) — (0(O)|T|0(O))^ — 



(19) 


( 20 ) 


5. An invariant in the probability distribntion |0j (f)p of states 


The time development of the initial state defined in equation 19 will be 


3 — iXot 




^Wl0(o)) = + 


72 


{ai\Vn-l) + a2\Vn-2) H - h am\Vn-m)) (21) 


This time dependent state can yield different time dependent probability distributions 
over the vertices of the graph according to different values for the amplitudes aj. Each 
Pj (t) = \(j)j (f)p will be either constant or periodic with a period of 7^. An invariant 
among these distributions will be the sum of the square of the difference between the 
maximum and minimum of each Pj [t), 


n—1 

^ ^ ^pmax _ pmin'^‘ 

3=0 


4 

n 


( 22 ) 


for any simply connected graph of n vertices. The proof is as follows; For any simply 
connected graph with n vertices the normalized eigenvector corresponding with the 
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eigenvalue Aq = 0 will be of the form 


l^^o) = 

/n 


/1\ 

1 


Vv 


The arbitrary unit vector in the eigenspace of \max can be represented as 


|aOp + joip + ■ ■ ■ + lUn-lP — 1- 


Equation can then be written 
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where Tq + + ■ ■ ■ + r‘^_^ = 1. 

Each (pj (t) will have its greatest or lowest norm when the corresponding 
term equals +1 or —1 respectively. 


+ 


1 Ti 


The maximum and minimum probability values and are then 

1 T ■ 

P— = |0+|2 = + ^ ^ 

■‘In \/n 1 


1 T ■ T- 
pmin _ 1^—|2 __j_ 3_ 

^ ' 2n ^ 2 


n—1 

j=0 


( Umax _ pmin\‘^ 
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j ) 
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( 23 ) 

( 24 ) 

( 26 ) 

.axt~\~^j ) 

(26a) 

(26b) 

(27a) 

(27b) 

(28) 


such that 




















CTQW amplitudes and invariants 

References 


10 


[1] S.E. Venegas-Andraca. Quantum walks: a comprehensive review. Quantum Information 

Processing, 11(5):1015-1106, 2012. 

[2] A. Ambainis. Quantum random walks, a new method for designing quantum algorithms. In 

V. Geffert and et al, editors, SOFSEM 2008: Theory and Practice of Computer Science, volume 
4910, pages 1-4. Springer Berlin/Heidelberg, Jan. 19-25 2008. 

[3] A. Ambainis. Quantum search algorithms. SIGACT News, 35(2):22-35, 2004. 

[4] A. Ambainis. Quantum walks and their algorithmic applications. International Journal of 

Quantum Information, 1:507-518, 2003. 

[5] E. Farhi and S. Gutmann. Quantum computation and decision trees. Phys. Rev. A, 58:915-928, 

1998. 

[6] H. Gerhardt and J. Watrous. Gontinuous-time quantum walks on the symmetric group. In 

S. Arora, K. Jansen, J. Rohm, and A. Sahai, editors, Approximation, Randomization, and 
Combinatorial Optimization.. Algorithms and Techniques, volume 2764 of Lecture Notes in 
Computer Science, pages 290-301. Springer Berlin Heidelberg, 2003. 

[7] A. Ahmadi, R. Belk, G. Tamon, and G. Wendler. Mixing in continuous quantum walks on graphs. 

Quantum Information and Computation, 3:611-618, 2003. 

[8] D. Ben-Avraham, E. Bollt, and G. Tamon. One-dimensional continuous-time quantum walks. 

Quantum Information Processing, 3:295-308, 2004. 

[9] A. Wyn-jones. Circulants. http://www.circulants.org/circ/, 2008. 

[10] P. Dukes. Quantum state revivals in quantum walks on cycles. Results in Physics, 4:189-197, 
2014. 



